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Abstract. The study of the content, distribution and kinematics of interstellar gas is a key 
to understand the origin and maintenance of both starburst and nuclear (AGN) activity in 
galaxies. The processes involved in AGN fueling encompass a wide range of scales, both spatial 
and temporal, which have to be studied. Probing the gas flow from the outer disk down 
to the central engine of an AGN host, requires the use of specific tracers of the interstellar 
medium adapted to follow the change of phase of the gas as a function of radius. Current 
mm-interferometers can provide a sharp view of the distribution and kinematics of molecular 
gas in the circumnuclear disks of galaxies through extensive CO line mapping. As such, CO 
maps are an essential tool to study AGN feeding mechanisms in the local universe. This is 
the scientific driver of the NUclei of GAlaxies (NUGA) survey, whose latest results are here 
reviewed. On the other hand, the use of specific molecular tracers of the dense gas phase can 
probe the feedback influence of activity on the chemistry and energy balance/redistribution 
in the interstellar medium of nearby galaxies. Millimeter interferometers are able to unveil 
the strong chemical differentiation present in the molecular gas disks of nearby starbursts and 
AGNs. Nearby active galaxies can be used as local templates to address the study of more 
distant galaxies where both star formation and AGN activity are deeply embedded. 



1. The feeding of nuclear activity in galaxies 

The phenomenon of nuclear activity in galaxies is viewed as the result of the feeding of 
supermassive black holes. Active Galactic Nuclei (AGN) are fed with material lying originally 
in the disk of the galaxy, therefore far away from the gravitational influence of the black 
hole. This implies in practice that the gas has to lose nearly all of its angular momentum 
on its way to the central engine. A correlation has been found between the presence of ~kpc- 
scale non-axisymmetric perturbations and the existence of activity in quasars, reckoned as high 
luminosity AGNs (HLAGNs). However, the search for a universal feeding mechanism in low 
luminosity AGNs (LL AGNs) had proved unsuccessful thus far [TJ El EJ 0] . The feeding problem 
is complicated by the fact that the AGN duty cycle (~10 6 ~ 7 years) might be shorter than 
the lifetime of the feeding mechanism itself (5j |6] . This time-scale conspiracy translates into a 
spatial-scale problem that can only be tackled by probing the critical scales for AGN feeding 



(< 10-100 pc). These spatial scales can only be reached with current mm- interferometers in 
nearby LLAGNs, a similar study in HLAGNs awaits the advent of future large interferometers 
like the Atacama Large Millimeter Array (ALMA). The gas masses involved in the feeding of 
HLAGNs, and possibly also the mechanisms of angular momentum transfer, could be different 
compared to those of LLAGNs [7J. A study of AGN feeding in nearby LLAGNs can help us 
to address the fueling problem in general, however. In particular, these type of studies lay 
the ground for future surveys where the whole luminosity range of AGNs will be explored. As 
molecular gas represents the dominant phase of neutral gas in galactic nuclei, rotational lines of 
carbon monoxide (CO) are well suited to undertake high-spatial resolution (— 1") interferometer 
mapping of the central kiloparsec disks of AGNs. 

1.1. The NUGA survey 

The NUclei of GAlaxies-NUGA-project, described by [BJ[9], is the first high-resolution (~0.5"- 
1") 12 CO survey of a sample of nearby (.D=4-40 Mpc) LLAGNs including the full sequence of 
AGN activity types (Seyferts, LINERs and transition objects). Observations have been carried 
out with the IRAM Plateau de Bure Interferometer (PdBI). The core sample of NUGA includes 
12 targets for which we have the highest spatial resolution (<10-100 pc) and sensitivity images 
(3-CT-detection limit~a few 1O 5_6 M ). This core sample is being expanded to reach a ~25 galaxy 
sample by including data of other LLAGNs from the PdBI archive. NUGA maps allow us to 
probe the gas flows at critical spatial scales where secondary modes embedded in the kpc-scale 
perturbations are expected to take over in the fueling process. To connect the large scales with 
the small scales, and have a complete picture of the fueling processes, a parallel HI survey has 
been conducted with the Very Large Array (VLA). This survey is used to probe the distribution 
and kinematics of neutral gas in the outer disks of NUGA galaxies [TUl [HI 12] . 

One of the main results of NUGA is the identification of a wide range of morphologies in 
the central kpc-disks of these LLAGNs. This includes one-arm spirals or m = 1 instabilities 
(NGC 4826: [DJ; NGC 3718: [13]), symmetric rings (NGC 7217: [H]; NGC 3147: [IS]), gas bars 
and two-arm spirals (NGC 4569: [17]; NGC 2782: [IE]; NGC 6574: [TO]; NGC 6951: [HI EI]), 
as well as mostly axisymmetric disks (NGC 5953: |16j). By itself, this result already suggests 
that several mechanisms, instead of one single universal mechanism, can cooperate to fuel the 
central engines of LLAGNs. 

1.2. Gravity torque maps: a new tool to quantify AGN fueling 

We have made a detailed case-by-case analysis of the distribution and kinematics of molecular gas 
in the galaxies of the core sample of NUGA, and interpreted these in terms of evidence of ongoing 
feeding. Our case-by-case in-depth study is the approach adopted to take advantage of the high 
quality of the NUGA mapall]. The gas response to the stellar potential is characterized with the 
help of high-resolution optical and NIR images of the galaxies. More precisely, we quantify the 
efficiency of the stellar potential to drain the gas angular momentum in a particular galaxy by 
deriving the gravity torques exerted by the potential on the gas [20] . With this method we can 
look for the feeding agent in the stellar potential. 

The results obtained from the analysis of stellar torques have revealed a puzzling feeding 
budget in the circumnuclear disks of the first NUGA targets studied by [2D]. Paradoxically, 
feeding due to the stellar potential seems to be presently inhibited close to the four AGNs 
analyzed, as torques are seen to be positive at r < 200 pc. To solve the puzzle, [2DJ suggest that 
gravity torques could be assisted by other mechanisms that become competitive in due time with 
non- axisymmetric perturbations. The authors estimate on a case-by-case basis that the gravity 
torque barrier associated with the ILRs of the galaxies analyzed in their paper could be overcome 



1 A statistical analysis of the whole ~25 galaxy sample is currently underway 




Figure 1. Snapshots of the stellar (left panel) and gas distributions (mid and right panels) 
at different epochs (in Myr) as seen in the numerical simulations that follow the evolution of 
a model galaxy to fit the distribution and kinematics of molecular gas in the NUGA galaxy 
NGC 2782 [18j . Axes are in kpc. Through dynamical decoupling of the nuclear bar with respect 
to the outer bar, the gas falls progressively inward. The gas infalling along the "dust lanes" is 
aligned along the nuclear bar, where the CO is observed in NGC 2782. Figure adapted from 



by viscosity (a source of negative torques), provided that the bar strength is moderate-to- low. 
Gravity torques and viscosity could thus combine their efforts to produce recurrent episodes 
of activity during the typical lifetime of any galaxy, a situation reminiscent of a self-regulation 
scenario, commonly invoked to understand the laws governing star formation and gravitational 
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Figure 2. a) (Upper panel) We overlay the 
HI contours with the the gravity torque 
map of NGC 4579. The derived torques 
change sign as expected if the butterfly 
diagram, defined by the orientation of 
quadrants I-to-IV, are due to the action of 
the large-scale bar. h)(Mid panel) The 
same as a) but for CO(l-O). c)(Lower 
panel) The same as a) but for CO(2-l). 
Figure adapted from [T2] . 
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Figure 3. The average fraction of the angular 
momentum transferred from/to the gas in one rotation 
(dL/L) are plotted as a function of radius, as derived 
from the HI (upper panel), CO(l-O) (mid panel) and 
CO(2-l) (lower panel) maps of the disk of NGC 4579. 
The locations of BAR resonances (ilLR, oILR, CRbar, 
VRRbar) and those of OVAL resonances (CHoval? 
UHRqvml) are highlighted. Figure adapted from 



instabilities in galaxy disks. The self-destruction of bars, triggered by the radial re-distribution 
of gas and the implied angular momentum exchanges, illustrated through numerical simulations 
|22l [23] , may have important side-effects for AGN fueling: several AGN feeding episodes may 
occur during a single bar cycle [20] . 

The existence of short-lived (lifetime < 10 6 -10 7 yr) feeding agents in the stellar potential 
(such as nuclear bars, oval perturbations or m = 1 instabilities), may also explain the above 
mentioned gravity torque puzzle. In particular, this solution would account for the low chances 
of finding a smoking gun evidence of fueling, especially if the number of targets scrutinized is 
small [2"U] . 

As recently shown |18[ 112) . by analyzing the cases of the NUGA galaxies NGC 2782 and 
NGC 4579, to find a galaxy caught in the act of fueling the AGN is unlikely but not impossible, 
however. The disk of NGC 2782, analyzed by [18], has two embedded bar-like perturbations: a 
primary large-scale oval (of 5 kpc diameter) and a prominent nuclear bar (of 1.5 kpc diameter). 
The decoupling of the nuclear bar has made it possible to transport molecular gas inside the ILR 
of the primary bar and then fuel the nuclear starburst and possibly the AGN. N-body simulations 
predict the decoupling of the nuclear bar and the ensuing gas inflow down to r ~200 pc (Fig. [I]). 
More compelling, in terms of AGN feeding evidence, is the case of NGC 4579, recently revisited 
in [12]. NGC 4579 has a prominent bar+spiral structure in the outer disk, r >2 kpc, as well 
as an oval distortion in the inner r ~200 pc. The gravity torque budget in NGC 4579, derived 
by [12] , shows that the fueling process is at work on different spatial scales (Figs. [2] and [3]). In 
the outer disk, the decoupling of the spiral structure produces net gas inflow on intermediate 
scales. Most remarkably, the co-rotation barrier seems to be overcome due to secular evolution 
processes, a result found in other NUGA galaxies [TT]. The gas in the inner disk of NGC 4579 
is being funneled down to r ~300 pc. Closer to the AGN (r <200 pc), gas feels negative torques 
due to the combined action of the large-scale bar and the inner oval. The two m = 2 modes act 
in concert to produce net gas inflow down to r ~50 pc, providing a clear smoking gun evidence 
of fueling with associated short dynamical time-scales (~l-3 rotation periods). 

The future instrumentation, of which the ALMA interferometer is a good example, will allow 
us to improve the statistics of these type of AGN surveys, by increasing the sizes of galaxy 
samples by orders of magnitude. In addition, the spatial resolutions typically reached by ALMA 
will provide a sharp view of the distribution and kinematics of molecular gas in the central pc 
regions of many nearby AGNs. Finally, we will also be able to analyze in detail the bar/AGN 
feeding cycles of distant galaxies and therefore compare fueling mechanisms as a function of 
redshift. 

2. The feedback of activity in galaxies 

The use of tracers specific to the dense molecular gas phase is a key to probe the feedback 
of activity on the interstellar medium of galaxies in the local universe. With typical densities 
n {H2) >10 4-5 cm -3 ), the dense molecular gas phase is the component directly involved in the 
fueling of star formation and AGN episodes. The onset of activity in its different varieties 
(starburst or AGN) can cause the injection of vast amounts of energy into the molecular gas 
reservoirs of galaxies through strong radiation fields and massive gas flows. The availability of 
high-spatial resolution observations of molecular species other than CO makes it possible to use 
the study of molecular gas chemistry as a tool to track down galaxy evolution. In this context, 
it is worth noting that extragalactic chemistry can put major constraints on the predictions of 
current chemical models of molecular clouds. For historical reasons, these models have been 
confronted thus far to observations of molecular clouds in our Galaxy. Active galaxies can drive 
chemical complexity to a higher degree and also on larger scales compared our Galaxy, however. 
More important, extragalactic chemistry studies are a key to constrain conversion factors, which 
are required to derive gas masses from line luminosities. This is mostly relevant considering that 



the different phenomena associated to activity are suspected to heavily influence the excitation 
and/or the chemistry of many of the tracers of the dense molecular gas phase which are routinely 
used in extragalactic research. This information is crucial for any quantitative assessment of the 
efficiency of the star formation process in galaxies. 

2.1. Molecular gas chemistry in starbursts 

Star formation can drive a large variety of chemical environments. Different processes can 
shape the evolution of molecular gas along the evolutionary track of a starburst: large-scale 
shocks, strong UV-fields, cosmic-rays, and X-rays. These processes are expected to be at work 
at different locations in the disk and/or the disk-halo interface of a starburst. During the 
fueling stage, gravity torques produced by density waves (bars and spirals) can funnel the gas 
inwards into the circumnuclear regions of galaxies. During this process, large-scale shocks can be 
generated in molecular gas. Furthermore, strong UV-fields are associated with the numerous HII 
regions of a typical starburst. The influence of UV fields in the surrounding medium can produce 
the propagation of Photon Dominated Region (PDR) chemistry in molecular gas disks. In a 
more advanced stage of a typical starburst episode, different types of chemistry can propagate 
into the halo of galaxies, entrained by massive galactic outflows. 

Different strategies can be adopted to tackle extragalactic chemistry. Single-dish line surveys 
allow the frequency coverage in the search of complex molecules to be expanded by significant 
factors and can thus provide us with a global (low-spatial resolution) view of the molecular 
gas inventory of a starburst. The single-dish line survey of the prototypical starburst galaxy 
NGC 253 is a good illustration of this technique [24"1 l25l 126] . 

The use of interferometers, with their high spatial resolution capabilities, makes possible 
to make a giant step forward in the field. The examples of M 82 and IC 342, two case study 
starbursts extensively observed with the PdBI and the Owens Valley Radio Observatory (OVRO) 
interferometers, are paradigmatic in this respect. Virtually all of the large-scale SiO emission 
detected in the PdBI map of M 82 [27] traces the disk-halo interface of the galaxy where episodes 
of mass injection are building up the gaseous halo. In contrast, widespread emission of the formyl 
radical, HCO, mapped in M 82 with the PdBI, reveals the propagation of PDR chemistry inside 
the disk of this starburst [28J. As such, the markedly different distributions of SiO and HCO 
in M 82 revealed by the PdBI, illustrates the power of mm-interferometry to unveil the strong 
chemical differentiation present in the molecular gas disks of starbursts. The scenario of a giant 
PDR in the disk of M 82 has received further observational support from the high abundances 
derived for molecular tracers which are specific to PDR environments such as CN, HOC + , CO + 
|29} 130] . Maps recently obtained in CN and HOC + lines with the PdBI are being used to 
complete this modelling effort |311 132j . In the case of IC 342, marked differences in morphology 
are also revealed among the different molecular maps obtained |33|, [34"] , using the OVRO and 
PdBI interferometers. Whereas the emission of tracers of PDR chemistry, like C2H and C 34 S, 
originates exclusively from the inner 50-100 pc circumnuclear region, where UV fields are high, a 
different group of molecules, including SiO, CH3OH, and HNCO, is correlated with the expected 
location of large-scale shocks induced by the bar potential in the outer disk. 

2.2. Molecular gas chemistry in AGNs 

Based on the observation of local template galaxies, it has been suggested that the excitation 
and/or the chemistry of molecular gas can show significant differences between starbursts and 
AGNs |35[ 136] . A recent study has shown that starbursts and AGNs populate different 
regions in a set of diagnostic diagrams that make use of several HCN and HCO + line ratios. In 
particular, Large Velocity Gradient (LVG) fits to the observed ratios indicate lower H2 densities, 
and larger HCN abundances in AGNs compared to starburst galaxies. Relative to other tracers 
of the dense molecular gas, overluminous HCN lines seem to be common in the circumnuclear 
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Figure 4. (Left panel) The CO(l-O) integrated intensity map of the Seyfert 2 galaxy NGC 1068 
|47] using the PdBI. (Right panel) Same for the SiO(2-l) line, but here showing the inner 
r ~300 pc region of the circumnuclear molecular disk (CND) of the galaxy mapped with the 
PdBI. The filled ellipse represents the beam size. Figure adapted from 
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Figure 5. (Left panel) The HCN(l-O) line intensity map obtained with the PdBI towards the 
CND of NGC 1068. (Middle) and (right panels) The HCO+/HCN, HNC/HCN and CN/HCN 
line ratio maps of the CND of NGC 1068 derived from PdBI data. Contours in the right panel 
stand for X-ray emission mapped by Chandra from 0.25 to 7.50 keV [50J. Figure adapted from 
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disks (CND) of many Seyfert galaxies (37J E3 EEl [36]. It is tempting to link the existence of 
overluminous HCN lines with the influence on molecular gas of some of the distinctive ingredients 
of AGN activity. Among these, X-rays are suspected to be able to process large column densities 
of molecular gas around AGNs, producing X-ray dominated regions (XDR). The abundances 
of certain ions, radicals and molecular species can be enhanced |39} 14*01 [4"T| I42j . In particular, 
and while still a highly debated issue, it has been discussed that the abundance of HCN can be 
boosted in XDRs, partly accounting for the reported observational dichotomy between starbursts 
and AGNs. Alternatively, the excitation of HCN lines in AGNs, some of them being sources 
with a high a infrared luminosity, might be affected by IR pumping through a 14^zm vibrational 



transition of the molecule |43^ HI] . 

The large HCN abundances measured in the nucleus of the Seyfert 2 galaxy NGC 1068 were 
the first observational evidence that molecular gas chemistry can be influenced by AGN activity 
|37|, 145] . The IRAM 30m telescope was used [38J to make complementary observations of the 
CND of NGC 1068 in eight molecular species. The global analysis of the survey, which includes 
several lines of SiO, CN, HCO, H 13 CO+, H 12 CO+, HOC+, HCN, CS, and CO, suggests that the 
bulk of the molecular gas emission in the CND of NGC 1068 can be interpreted as coming from a 
giant XDR. Recently, a high-resolution PdBI map has provided a sharp view of the distribution 
of the SiO emission coming from the CND of NGC 1068 (Fig. The enhancement of 
SiO in the CND of NGC 1068 cannot be attributed to the action of ongoing star formation as 
there is negative evidence of a recent starburst in the circumnuclear region of this Seyfert [15] . 
Alternatively, the processing of 10 A dust grains by X-rays, as a mechanism to enhance silicon 
chemistry in gas phase, might explain the large SiO abundances of the CND. The inclusion 
of dust grain chemistry could solve the controversy between different gas-phase XDR schemes 
regarding the abundance of some molecular tracers like HCN in X-ray irradiated environments. 
New maps of NGC 1068 obtained with the PdBI in molecular species like HCN, HCO+, HNC 
and CN are being used to study the chemical differentiation inside the CND of the galaxy, and 
explore its dependence with the illumination of molecular gas by the strong X-ray source of the 
nucleus 09] (Fig. [5}. 

2.3. Molecular gas chemistry in IR luminous galaxies 

The use of HCN lines as a quantitative tracer of the dense molecular gas in galaxies can be 
questioned if the excitation conditions and/or the chemical environment of molecular gas depart 
from normality. This can alter the conversion factor between the luminosity of HCN lines and 
the mass of dense molecular gas. The doubts cast on HCN call for the observation of alternative 
tracers of dense gas in luminous and ultraluminous infrared (IR) galaxies (LIRGs and ULIRGs). 
This question is crucial to disentangle the power sources of the infrared luminosities of these 
galaxies, where both star formation and AGN activity are embedded. 

A recent paper [5T] has presented evidence that the Lfir/Lhcn(i-o) ratio, taken as a proxy for 
the star formation efficiency of the dense gas (SFEdense)) is a factor ~2-3 higher in IR luminous 
galaxies (-Lfir > 10 11 L Q ) compared to normal galaxies. Local universe LIRGs and ULIRGs 
populate a region in the SFEd en se diagram that lies between those occupied by normal and high-z 
IR luminous galaxies. The reported trend in the SFEd en se derived from HCN data implies that 
there is a statistically significant turn upward, in the KGnnicutt-Schmidt law, 5]sfr ^ dense' ^ 
high L Fm : N changes from ~0.80-0.95 (for L Fm < 10 11 L Q ) to ~1.1-1.2 (for L Fm > 10 11 L Q ). 
These results confirm the predictions of star formation models [52]. In addition, and based on a 
multi-line LVG analysis of HCN and HCO + data that follows the survey published in [53], [51] 
find that the the conversion factor between -^hcn(i-o) an d the mass of dense molecular gas is 
~3 times lower at high Lfir- Both factors taken into account, this result reinforces a scenario 
where the SFEd cn se could well be up to an order of magnitude higher in extreme LIRGs/ULIRGs 
compared to normal galaxies. Of particular note, a significant overabundance of HCN has also 
been reported in the IR luminous z~4 quasar APM 08279 [54, 44J, an example of extreme 
chemistry and excitation conditions of molecular gas at high redshift. 
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